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Abstract 

We calculate the dissociation rates of quarkonium ground states by tunnelling and 
direct thermal activation to the continuum. For hadronic matter at temperatures T < 
0.2 GeV, neither of these mechanisms leads to a sufficiently large dissociation to explain 
the experimentally observed suppression of charmonium. Dissociation by sequential 
excitation to excited energy levels, although OZI-forbidden, requires further analysis. 
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For sufficiently heavy quarks, the dissociation of quarkonium states by interaction 
with hght hadrons can be fully accounted for by short-distance QCD [1-4]. Pertur- 
bativc calculations become valid when the space and time scales associated with the 
quarkonium state, rg and tq, are small in comparison to the nonperturbative scale 

« ^QCD' (la) 

tQ « Aq^d; (1^) 

Aqcd ^^^^ characteristic size of the light hadrons. In the heavy quark limit, the 
quarkonium binding becomes Coulombic, and the spatial size tq ~ (agmg)"^ thus is 
small. The time scale is by the uncertainty relation given as the inverse of the binding 
energy Eq ~ mg and hence also small. For the charmonium ground state J/i/j, we have 

~ 0.2 fm = (1 GeV)-^ (2) 



and 

= 2Md - ~ 0.64 GeV. (3) 

With Aqcd — 0.2 GeV, the inequalities (1) seem already reasonably well satisfied, and 
also the heavy quark relation E^ — (l/mgr^) is very well fulfilled. We therefore expect 
that the dissociation of J/ip's in hadronic matter will be governed by the J/V'-hadron 
break-up cross section as calculated in short-distance QCD. 

Nevertheless, in view of the finite charm quark mass, it makes sense to study 
possible non-perturbative contributions, in particular for the calculation of the break- 
up process in matter, where in addition thermal activation may come into play. For an 
isolated J/i/^-hadron system, non-perturbative interactions can be pictured most simply 
as a quark rearrangment. Consider putting a, J /if; "into" a stationary light hadron; the 
quarks could then just rearrange their binding pattern to give rise to transitions such as 
J/ip + N ^ Ac + D OT J/ip + p ^ D + D (Fig. 1). The probability for such a transition 
can be written as 

Prearr- J dhR{r) 10^.(01', (4) 

where the spatial distribution of the cc bound state is given by the squared wave function 
|(/)^(r)p. The function R{r) in cq. (4) describes the resolution capability of the colour 
field inside the light hadron. Its wave length is of order Aqqj-j, and so it cannot resolve 
the charge content of very much smaller bound states; in other words, it does not "see" 
the heavy quarks in a bound state of radius tq « Aq^j-, and hence cannot rearrange 
bonds. The resolution R{r) will approach unity for tAqcd >> 1 and drop very rapidly 
with r for tAqcd < 1, in the functional form 



i?(r) ~ (rAQCD)^ rAQCD<l, (5) 

with n = 2 [5] or 3 [1]. As a result, the integrand of eq. (4) will peak at some distance 
ro, with Tq < Tq < Aqcd- Since the bound state radius of the quarkonium ground state 
decreases with increasing heavy quark mass, while R{r) is mg-independent, — > rg — > 



1 



as rriQ — > oo. Hence Pr vanishes in the hmit mq — > oo because R{ro) does, indicating 
that the hght quarks can no longer resolve the small heavy quark bound state. 

In a potential picture, the situation just described means that the charm quarks 
inside the J/i/) have to tunnel from r = r^p out to a distance at which the light quarks 
can resolve them, i.e., out to some r ~ cAqqj^, where c is a constant of order unity (Fig. 
2). Such tunneling processes are therefore truly non-perturbative: they cover a large 
space-time region, of linear size ^qcd^ involve any hard interactions. The 

first aim of this letter is to estimate the contribution of non-perturbative tunnelling 
to the dissociation of quarkonium states. Following this, we turn to matter at finite 
temperature, where the J/ijj can be excited into the continuum by thermal activation. 
Our second objective is to calculate the rate of dissociation by this mechanism. 

In general, the problem of quark tunneling cannot be solved in a rigorous way, 
since it involves genuine non-perturbative QCD dynamics. However, the large mass of 
the heavy quark allows a very important simplification, the use of the quasiclassical 
approximation. In this approximation, the rate of tunneling i^tun can be written down 
in a particularly transparent way: it is simply the product of the frequency a;^ of the 
heavy quark motion in the potential well and the tunnelling probability Ptun when the 
quark hits the wall of the well, 

(6) 

The frequency a;^ is determined by the gap to the first radial excitation, 

oj^ ~ (M^/ - M^) ~ E^. (7) 

Consider now the potential seen by the cc (Fig. 2). For a particle of energy E, the 
probability of tunneling through the potential barrier V{r) is obtained from the squared 
wave function in the "forbidden" region. It can be expressed in terms of the action W 
calculated along the quasiclassical trajectory, 

Ptun = e-2^, (8) 

where 

W= \p\ dr, (9) 
J n 

Here the momentum \p\ is given by 

\p\ = [2MiVix)-E)]'/\ (10) 
and ri , r2 are the turning points of the classical motion determined from the condition 

vin) = E. 

In our case, the width of the barrier is approximately 0.6 Aq^j^, while its height 
iy — E) is equal to the dissociation threshold Eq. The mass M in Eq.(lO) is the reduced 
mass, M = mQ/2. We thus have 

W ~ 0.6 y/ruQ Eq/Aqcd. (11) 
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For the J/V', we get from (11) the value ~ 3; this a posteriori justifies the use of 
quasiclassical approximation, which requires S > 1. 

Using eq. (11), we obtain as final form for the tunneling rate (6) 



Rtun = exp - (1.2^mc E^/Aqcd)- (12) 

With the above mentioned J/ij; parameters, this leads to the very small dissociation 
rate 

Rtun ^ 9.0 X 10-3 fm-^ (13) 
In terms of Rq, the J/iIj survival probability is given by 

-^tun = exp <^ - / dt Rtun >, (14) 



where tmax denotes the maximum time the J/tl) spends adjacent to the light hadron. 
In the limit tmax oo, S^p vanishes. However, the uncertainty relations prevent a 
localisation of the two systems in the same spatial area for long times. From Ax < Aq^q 
we get Ap > Aqcd, so that the longest time which the J ftp can spend in the interaction 
range of the light hadron is 



2 N^/' 
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with m for the mass of the light hadron. For nucleons or vector mesons, this time is 4 
- 5 fm, and with this, the survival probability is very close to unity, so that tunnelling 
cannot result in a noticeable non-pcrturbativc J/ijj dissociation. 

In a medium at finite temperature, however, the J/i(j can in addition be thermally 
excited into the continuum and thus become dissociated; this thermal activation could 
be non-perturbative. Here we shall simply consider a J/V^ in a thermal medium of 
temperature T and calculate its excitation rate, without asking how the constituents of 
this medium bring the ground state cc into the continuum. Since hadrons in a medium 
of temperature T may not be able to interact sufficiently hard with the J/ i/j to overcome 
the mass gap to the continuum, we obtain in this way an upper bound to direct thermal 
dissociation. 

The partition function of the system at finite temperature is 

Z(T) = e'^"'''^ = ^^(^) + ^cont(T), (16) 
n 

where .^cont and are the continuum and the boimd-state contributions, respectively. 
We assume that the cc will be distributed among the ground state J/i(^ and the contin- 
uum above E = 2Md according to 

piE,T) = c{E) e-^/^[5(l - E/M^) + Q{E - 2Md)], (17) 



where c{E) denotes the degeneracy factor, which we take to be constant, c{E) = c. The 
continuum part of the partition function is given by 

ZcontiT) = Vcj ^e-^/^0(E - 2Md), (18) 

where V is the canonical volume containing the system. Changing variables to 

p^dp = MV2MEdE, (19) 



we get 

^cont(T) = Vcj^,M ^ T^l^ e-^^l\ (20) 

while the bound-state part is obtained from eq. (17) as 

Z^(T)=c, (21) 

with E^ = 2Md - and c = c e'^^/^. 

The rate of escape into the continuum can be estimated by the average time needed 
to leave the spatial region of the potential well, 

R,ct^{t)-'^{viE))/L. (22) 

Here v{E) is the velocity of the (reduced) charm quark in the continuum, and L ~ 
(1 — t^Aqcd) -Aqcd distance from the average J/ip radius to the top to the 

potential well. Changing variables from the overall energy to the energy above and 
using v{E) — p{E)/M, we bring eq. (22) into the form 

_ 1 V f d^p 

-n-act — 



Z(T)MLj (27r)3 



/ 7|;;f3Wce(E-^^)e-^/^]. (23) 



With the substitution p{E) — ^J2M{E — E^) we then have after integration the result 

1 V 



Ract — 



Z{T) L 



{^MT^^ e-^*/^, (24) 



where V — is the canonical volume. 

For high temperatures, T ^ E^, the continuum gives the dominant contribution 
to the statistical sum, so that we can replace Z{T) by Z^ont in eq. (24) to get 



Recalling that the thermal velocity of a free particle in three dimensions is just Vth{T) = 
we thus recover the classical high-temperature limit for the thermal acti- 
vation rate 

B.C. = (26) 
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At low temperatures, for T <^ E^^ the discrete term in p give the main contribution to 
Z{T). In this temperature range, which is the one of interest here, we thus obtain from 
eqs. (24) and (21) 

i?act = ^ M e-^^/^, (27) 

as our final result. 

The explicit appearance of the volume factors L in the above results for thermal 
activation may at first sight seem strange. Its origin is the fact that although the 
intrinsic Boltzman factor for excitations to the continuum is small, the density of states 
there is quite large. Nevertheless, the limit of large L is well defined, since for L ^ oo 
the continuum part (20) of the partition function dominates, so that the result (27) 
is replaced by the classical formula (25). The heavy quark limit M ^ oo is also well 
defined, since at very large M the size of the system shrinks to L ~ 1/M. It is moreover 
important to remember that for heavy quark- antiquark systems the binding energy (on 
which the rate (27) depends exponentially) increases with the mass of the quark. 

With the same charmonium parameters as above (£/^ — 0.64 GeV, L — 0.6 Aqqj-,, 
2M = rric = 1.5 GeV) and at the temperature T = 0.2 GeV, we thus obtain 

i?act - 6 X 10-3 fm-\ (28) 

which is of the same size as the tunnelling rate (13). The precise value of L is somewhat 
uncertain, of course. We feel, however, that L ~ 1 fm is an upper bound, leading to 
i?act < 2 X 10-2 fm-^ as upper bound for the rate. Thus direct thermal activation can 
also not provide a significant amount of J/ifj dissociation at temperatures below the 
binding energy. Already for static media of 4-5 fm life-time the survival probability is 
very close to unity; any expansion would further reduce activation effects. 

We should however be cautious about the application of the thermal activation for- 
mula used above. Because of the large density of states in the continuum, the treatment 
of the factor L was somewhat heuristic. To better understand this dependence, it will 
probably be necessary to investigate directly the sequential transitions from the ground 
state to excited states and then to the continuum. Thus the J/^/J could be dissociated 
by first exciting it to a Xc and then bringing the Xc to the continuum. Since the phase 
space for the continuum is so large, the probability of de-excitation back to a bound 
state is negligible. Unfortunately, it is difficult to obtain a simple estimate of the various 
processes which contribute to this sequential excitation. The analysis of this process 
awaits further work. 
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Figure Captions 



Fig. 1 Rearrangement transition p + J/\E' — > D + 
Fig. 2 J/^f dissociation by tunneling. 
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